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In loop quantum cosmology, the Hamiltonian reduces to a finite difference operator. We study the 

initial singularity and the large volume limit against the ambiguities in the discretisation and the 

operator ordering within a homogeneous, isotropic and spatially flat model with the cosmological 

constant. We find that the absence of the singularity strongly depends on the choice of the operator 

ordering and the requirement for the absence singles out a very small class of orderings. Moreover 

J— ^ ' we find a general ordering rule required for the absence of the singularity. We also find that the 

*vj , large volume limit naturally recovers a smooth wave function in the discretisation where each step 

■ corresponds to a fixed volume increment but not in the one where each step corresponds to a fixed 

i^^' area increment. If loop quantum cosmology is to be a phenomenological realisation of full loop 

quantum gravity, these results are important to fix the theoretical ambiguities. 

PACS numbers: 04.60. Kz, 04.60. Pp, 98.80.Qc 
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I. INTRODUCTION 



{3J[)| Loop quantum gravity (LQG) is a background-independent nonpcrturbative quantum gravity. There, one quantises 

the Hamiltonian formulation of general relativity based on the SU(2) connection A^^ and the densitised triad Ef on the 
three-dimensional space, where a, 6, ... and i, j, ... both run over 1, 2 and 3. a, 6, ... are tangent indices, while i, j, ... 
" ►^ ■ are indices associated with a basis {ri} for the Lie algebra of SU(2) Lie group and raised and lowered with Kronecker's 
QQ ■ delta Sij. The three-dimensional inverse metric is given by q"-'' = \det E'^\~^E°-E''S'^^ . To quantise the system, LQG 

^~~' invokes holonomies he = eJ='*<»'^*''^ along an edge e and fluxes Eg — J^cPanaEfr^ over a two-surface S, where (Pa 

jT •, , is the area element on S and Ua is a unit normal to S. There appears a nondimensional constant parameter 7, 

. ' which is called the Barbero-Immirzi parameter and cannot be determined within the theory. The kinematical Hilbert 

1^ ' space is spanned by the spin network states. The Hamiltonian as well as the area and the volume are constructed 

from holonomies and fluxes and act as operators on the kinematical Hilbert space. The Hilbert space is spanned by 

equivalence classes of the spin networks under diffeomorphisms or the s-knot states. See [i|-y] for the basics and recent 

developments in LQG. 

Loop quantum cosmology (LQC) [^ is a quantised minisuperspace model motivated by LQG. In traditional quantum 
cosmology [5j, symmetry reduced models are quantised in the usual Schrodinger representation and the quantised 
5^ Hamiltonian constraint yields the Wheeler-DeWitt (WDW) differential equation, while in LQC one uses the so-called 

Ci ' polymer particle representation [^ which is unitary inequivalent to the Schrodinger representation and obtains a 
second order difference equation rather than the differential equation. Remarkably, Bojowald T\ demonstrated that 
there can be no big bang singularity in the following two aspects. First, the spectrum of the inverse scale factor 
operator is bounded from above. Second, the wave function of the universe can be uniquely extended beyond the 
point which was the initial singularity in classical theory. The above features have been first shown for a homogeneous, 
isotropic and spatially flat model and subsequently generalised to nonflat or anisotropic cases [4], although it is yet 
uncertain whether this can be generalised to inhomogeneous models (see e.g. [1, Q)- 

In this paper, we visit the absence of the singularity in the presence of the cosmological constant against the 
quantisation ambiguities. In particular, we focus on the choice of the discretisation parameter A which corresponds to 
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the coordinate length of edges for the basic holonomies and the operator ordering in the Hamiltonian constraint. To 
study these issues, we fix the cosmological model to be homogeneous, isotropic and spatially flat with the cosmological 
constant. We first adopt the constant A as in [3], for which the universe gains a fixed area quantum at each step and 
find that the absence of the singularity strongly depends on the operator ordering in the Hamiltonian. However, for 
this discretisation, we encounter a serious problem in the physical interpretation of the obtained wave function in the 
large volume limit. The problem is that the behaviour of the wave function obtained from the difference equation 
in LQC does not agree with that obtained from the WDW equation in the large volume limit. Since the effects of 
discreteness are dominant in the Planck scale physics but should disappear when the universe becomes large, the 
discrete wave function in LQC in the large volume is expected to have a correspondence to the smooth wave function 
which is a solution to the WDW equation. We see that this problem is resolved if A is chosen to vary [l^l so that the 
universe gains a fixed volume quantum and find that the absence of the singularity strongly depends on the operator 
ordering in the Hamiltonian also in this discretisation. It is already studied in detail how the large volume limit 
depends on the choice of the discretisation in the presence of matter fields by Nelson and Sakellariadou [111, [3 and 
the present result with the cosmological constant but without matter fields is consistent with theirs. 

This paper is organised as follows. In Section [Hi we introduce the discretised Hamiltonian for a homogeneous, 
isotropic and spatially flat model with the cosmological constant. In Section IIIIl we choose A to be constant and 
derive the difference equation. Then, we demonstrate how the absence of the singularity depends on the choice of 
the operator ordering and numerically and analytically show the large volume limit is problematic. In Section IIVI 
we vary A appropriately and derive the difference equation. Then, we show that the large volume limit problem is 
resolved in this discretisation and also find that the absence of the singularity strongly depends on the choice of the 
operator ordering. Section |V] is devoted to discussion for the influence of the matter fields. In Section FVll we conclude 
the paper. In this paper we use the units in which c = h ~ 1. 

II. LOOP QUANTUM COSMOLOGY 

A. Hamiltonian constraint 

In the present paper, we focus on a homogeneous, isotropic and spatially flat universe. In classical theory, the line 
element for such a spacetime is given by the flat Friedmann-Robertson- Walker (FRW) metric 

ds^ = -dt^ + a{tf {dx^ + dy"^ + dz^) , (2.1) 

where a{t) is called the scale factor. To remove the divergence coming from the volume integral we introduce an 
elementary cell U on the three dimensional space. In this case, the gravitational Hamiltonian constraint is written 
as [11 

^- - 16^ / ''^ (-^^^^^^-^-^"^^" + ^1^^^^^\^) , (2.2) 

where G, A, eijk and F^^ = daA\ — dbA\ + e*j,A^A^' are the gravitational constant, the cosmological constant, the 
Levi-Civita symbol and the curvature associated with the connection A^, respectively. Although the cosmological 
constant might be an emergent object from some unknown effects, we incorporate it into the Hamiltonian as usual for 
simplicity. In the elementary cell V , we define a time-independent fiducial flat metric qab, an associated ortho normal 
triad {e°} and cotriad {ujI}. In the flat FRW universe, the connection ^^ and the densitised triad Ef are given by 

^1 = cV^'^'^'^a, Et=pV-^^'^^^qe1, (2.3) 

where q is the determinant of qabi c = V^ sgn{p)jda/dt, \p\ — Vq c? and Vq = / d^x^fq. Vq is thus the volume 
of the elementary cell with respect to the fiducial metric. In other words, \p\ is proportional to the physical area 
of the elementary cell and c is the conjugate momentum of p. The Poisson bracket between c and p takes the form 
{c,p} = 87rG7/3. Now, the holonomy ft,^ along an edge parallel to e° is given by 

;j(A)^gAcr,^ (2.4) 

where A is the coordinate length of the edge and {ji] are a basis for the Lie algebra of SU(2) Lie group. The flux Es 
is simply given by Es = pV^ As, where As is the area of the surface S [I3. Note that in defining the Hamiltonian 



constraint operator, one traces over SU(2) valued holonomies, and then there appears an ambiguity in choosing the 
irreducible representation to perform the trace. Here, we choose the spin J = 1/2 representation for simplicity. Then 
the holonomy becomes 

/if) = cos (Ac/2) I + 2 sin (Ac/2) n, (2.5) 

where I is the unit 2x2 matrix and {tj} are related to the Pauli matrices {<Ji\ via 2irj — Ui. 

We rewrite the Hamiltonian constraint in terms of the flux p and the holonomy h\, as in the full theory. For the 
triad part of the first term in the constraint (12. 2p , we obtain 
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where e"'"^ is the Levi-Civita symbol, V is the volume and {•, •} denotes the Poisson bracket. For the curvature F*j,, 



we rely on a standard prescription in gauge theories. We consider a loop a, which is a square Dy spanned by two 
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triad vectors ef and e^ of which each side is as long as A in the coordinate length. Then, the ab component of the 
curvature is given by 
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where the holonomy ha ij along a = D^ is the product of holonomies along the four edges. 






Substituting Eqs. (gH), ([221), and 
1 / 4sgn(p) 



into Eq. (|2.2p . Cgrav can be expressed as 
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where we have used the relation 
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To quantise the Hamiltonian, we replace h] and p with the corresponding operators. The kinematical Hilbert 
space is defined by H^^^^ — £^ (Rsohr , dMBohr ) 1 where MBohr is the Bohr compactification of the real line and d//Bohr 
is the Haar measure on it [i|, y, |a, [l3] • An orthonormal basis for the kinematical Hilbert space is given by a set of 



eigenstates {|/i)} of p, which satisfy the orthonormality relations (/xil/12) 
on the state \p) is given by 
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The action of the triad operator p 
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where ^pi := vG is the Planck length. That is, the eigenvalues of p are labelled by the dimensionless parameter /j,. 
The states {|/x)} are also eigenstates of the volume operator V = \p\^/'^: 
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Using Eqs. (12. 4p and (12. 9p and replacing the Poisson bracket with a commutator, after some calculation, we find 
1 
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where the operator ordering is fixed for simplicity. The ambiguity in the ordering will be discussed later. Note that A 
itself is an operator in general. It should be noted that in LQC we do not take the limit A ^- 0. Later we will discuss 
the physical motivation for this setting. 



B. Discretisation ambiguity 

In full LQG, the geometry is quantised through the area and volume operators. To define the area operator A(5) 
for a two-surface S, we partition S into N small two-surfaces {5„} such that U„iS„ ~ S. For sufficiently large N, 
there is {5„} such that no iS„ will contain more than one intersection with the graph F of the spin network |S). The 
sum over n reduces to a sum over the intersection points between S and F and is independent of N for sufficiently 
large N. Then, the action of the area operator becomes 

A(5)|S) = 4^7^li E \/2j"(j;" + 1) + 2jf Of + 1) - JIUI + 1)|S), (2.15) 

i 

where {i} label the intersection points between the graph F and the two-surface S, the indices u, d and t stand 
for the edges upward, downward and tangential to the S, respectively, and the positive half integers j", jf and j* 
are the spins of the links labelled by i [I^. Thus, there appears the smallest area A = 2^/3TT^lp^. For the volume 
operator V for a three-volume TZ, we take a similar strategy. We partition the three-volume 7?. into cubes {TZn} of 
the coordinate volume e^ and for sufficiently small e no cube will contain more than one node. The volume operator 
has a nontrivial action only on nodes and hence it will no longer depend on the value of e. Then, it turns out that 
the spin network state is an eigenstate of the volume operator and the eigenvalue is given by the sum over the nodes 
which are contained in the three-volume TZ and at least quadrivalent. Similarly, we define the quantised Hamiltonian, 
which has a nontrivial action only on nodes. Thus, for sufficiently small e, the action of the quantised Hamiltonian 
will not depend on e. 

In LQC, however, we leave the parameter A of the discretisation nonzero finite. In fact, we are forced to do so in 
the present formulation of LQC because there is no operator corresponding to c. On the other hand, the physical 
results seem to depend on the choice of A. From this point of view, there is no first principle within the formulation 
of LQC about how to choose the nonzero finite value for the discretisation parameter A. We should probably fix the 
discretisation parameter A (and other quantisation ambiguities) in LQC so that LQC can reproduce the features that 
full LQG should have. 

To see how A is kept nonzero, we see the following relation: 

\^\h['^ = ^lll,\X\h['\ (2.16) 

Note that \p\ is the physical area of the elementary cell. Since the curvature in the Hamiltonian invokes the holonomy 

ha ij along the square Oij, it would be reasonable to assume that each edge of this square is quantised so that the 
holonomy along each edge intersects the smallest area A of full LQG. This argument motivates us to choose A to be 
a constant /iq of order unity. This choice is adopted in Refs. [7|, Il5l . Il7| . The value for fi^ to fulfill this requirement 

exactly is 3v3/2. The basic operator which appears in the Hamiltonian is exp {ifioc/2) and this acts on |/^) as follows: 

e^^Waj^) ^l^ + ^o). (2.17) 

This means that the eigenvalues of /i for the states which appear in the Hamiltonian constraint are spaced at constant 

intervals. Note that as seen in Eq. (I2.17p . the operator exp (i/ioc/2) drags the state |/i) by the parameter length /ip 
along the vector d/d/i. Then, we can rewrite Eq. (|2.17p as 

evW2|^) =e''«(''/''^V)- (2-18) 

We call this choice of the parameter A the equi-area discretisation. 

However, this is not the only possible choice. Since the curvature in the Hamiltonian invokes the holonomy ha ij 
along the square D^ , it would also be reasonable to assume that the area of this square is set to be the smallest area 
A of full LQG. This motivates us to choose A to be a function A = /i(p) such that 

fl^\p\^A. (2.19) 

This gives another choice of the parameter A. In this case, from Eq. (|2.1ip we find that p, depends on /i as follows: 

M^H = ^. (2.20) 



Similarly to Eq. (|2.18p , we can rewrite the operator exp {ific/2) as 

gi;W2|/,) = e^^'^/'^^'V)- (2-21) 

If we introduce a variable v satisfying 



dv = -dn, (2.22) 



we can rewrite the right hand side of Eq. (|2.2ip as 



e''^\n)^e^\fi). (2.23) 



This means that if we adopt the u-representation of the wave function, the action of the Hamiltonian becomes simple. 
The explicit integration of Eq. (|2.22p using Eq. (|2.20p yields 

V = sgn{^l)K\^l\i , (2.24) 



where K = 2-\/2/(3\/3'\/3). We should note that v is proportional to the volume V. This choice is adopted in 
Ref. [lOl- In this case, it is more convenient to use the eigenstates \v) of the volume operator V to see the action of 
the Hamiltonian constraint. The eigenstate \v) satisfies 

V\v)^V,\v), (2.25) 

where 

K=(^7/^.)'''M. (2.26) 

The basic operator which appears in the Hamiltonian is e*'^'^/^ and its action on \v) is given by 

7f^\v) ^ \v + 1) . (2.27) 

Consequently, the Hamiltonian constraint involves the volume eigenstates of which the eigenvalues are equally spaced. 
In other words, the spectrum of the volume operator is distributed at equidistant intervals in volume. Thus, we will 
call this discretisation the equi- volume discretisation. Note that since A — p.{p) depends on p, A should be treated as 
an operator in quantum theory. 

III. EQUI-AREA DISCRETISATION 

A. Absence of singularity and operator ordering 

First we concentrate on the equi- area discretisation. The holonomy operator acts on |/i) as 

4''° V) = ^(Im + Mo) + Im - Mo))I + ^(l/i + /^o) - l/i - /^o))rfe. (3.1) 

Even within this discretisation, there are many possible choices of the operator ordering. To make the notation simple, 
we put 

9 . 2M0C~^ /XqC --;-- . paC /ioC ~ Moc,. . PoC ,„ „. 

i* = sm — — cos"^— — , l^£j = sm——Vcos—- cos— — Ksm— — , (3.2) 

2 2' 22 22' ^ ^ 

where F and EE relate to the squared holonomy haij (|2.8|) and the Poisson bracket hi{{hi)~^,V} in Eq. (|2.6|) . 
respectively. For simplicity, hereafter we fix the ordering for the contents of F and EE, and then there are two choices 
of the operator ordering in the Hamiltonian constraint operator as FEE or EEF. We first consider the operator 
ordering adopted in Ref. |7| : 

1 / 96»(sgn(p) ) 9^ , , 2 



Csrav = ,^ ,2 , —^^^FEE + 2-f' AV . (3.3) 



Then the action of the operator on a state |^) leads to the difference equation, i.e., (/i|Cgrav|^) = yields 



(3.4) 

where ^(/i) = (/i|^). If we interpret the triad coefficient p as an internal time, we can regard the difference equa- 
tion p. 41) as an evolution equation with respect to the discrete time. We can see that ^(0) disappears for /i = and 
±4/xo. For this reason, the solution can be uniquely extended beyond the classical singularity ^ — 0. That is, given 
some two initial data *(e + ANfio) and ^(e + 4{N + I)/xo) for e € (0, 4/zo) and a natural number N, one can uniquely 
determine the values ^(e + 4n/io) for n = 0, ±1, ±2, • • • . For the e = case, given some two initial data '^{ANiiq) and 
^(4(iV + l)^o), the difference equation (|3.4I) generates \I'(4n/zo) for n = 1, 2, • • • but the set of \E'(8/io) and *(4^o) 
does not generally satisfy Eq. p.4p . This means that ^(4A^/io) and *(4(iV + l)/zo) are constrained, and once this is 
satisfied, we can uniquely determine the values ^{in^o) for n = ±1, ±2, • • • . Thus, we can conclude that the system 
has no initial singularity in this operator ordering. Although ^(0) is left undetermined in Eq. (3.4), we can determine 
unambiguously the discrete evolution beyond fj, — 0. That is, the unique quantum evolution is not affected in spite of 
the undetermined value at // = 0, which was the initial singularity in classical theory. This ambiguity is fixed so that 
^(0) = in Ref. 7] and we also use this value in our numerical analysis. 

Here we should note the superselection of the kinematical Hilbert space. As seen from Eq. (|3.4p . a subspace He 
with a fixed parameter e S [0, 4/io) spanned by a basis {\fj,) : fi — e + 4n/io, n e Z} of the kinematical Hilbert space 
^kin^ are closed with respect to the action of the Hamiltonian constraint. More precisely, the kinematical Hilbert 
space is naturally decomposed into the sectors with respect to the action of the Hamiltonian constraint as 

-T/giav _ /^V\ nj 
"•kin — VC <^' 

te[QA^la) 

We refer to this decomposition as the superselection and to these subspaces as the superselection sectors. The 
superselection in LQC indicates that for each sector we can fix the value of e and work within this sector T-L^- Since 
there is no initial singularity for the sectors with e 7^ in this case, we only have to choose the sector with e = to 
study the presence or absence of the initial singularity at /x = 0. On the other hand, when we consider the large /i 
limit, the value of e does not affect the qualitative behaviour of ^(/i) and hence we can focus on the sector with e = 
again. This is also the case in the equi- volume discretisation. 
Next we consider the following operator ordering: 

1 /96i(sgn(p))-j^^" 2, 



Then, the difference equation is given by 

\-V,+,, + V,^,, I {^{11 + 4/.o) - 2vl/(^) + ^{^, - 4mo)) + ^-^^^^^^KV,■^{^l) = 0. (3.6) 

We can see that Eq. p. 61) becomes trivial for /i = and hence we cannot determine ^(— 4/io) from 4'(4/io) and ^(0). 
This fact indicates that we cannot determine all ^(/i) through /z = from the data ^(/i) for /i > 0. In this sense, 
the model contains the initial singularity at /.t = 0, beyond which the evolution cannot be uniquely extended. This 
means that the absence of the singularity depends on the choice of the operator ordering in the quantisation of the 
Hamiltonian. 

B. Large volume limit problem 

In this section, we shall see the large- /Lt behaviour of the wave function *(/i) determined by the difference equation 
and discuss its physical implication. Since the effects of the discreteness are dominant in the Planck scale physics but 
should disappear when the universe becomes large, in LQC the discrete wave function in the large volume is expected 
to be well approximated by a smooth wave function. If this naive expectation is valid, it is natural to think that the 
smooth wave function is described by a solution to the WDW equation which is obtained by quantising the system 
in the usual Schrodinger representation. We give a brief derivation of the WDW equation in Appendix [X] 

To see the large- /i behaviour of the discretised wave function ^(/i), we numerically solve the difference equation 
p.4p . We first choose the initial values ^(— 4/io) and ^(0), and then determine ^(4n/^o) for n = 1,2,... by the 



difference equation. Here we define the dimensionless cosmological constant as A = (167r/3)7'^Z|,[A. In the numerical 
calculation below, we set A/Xg = 0.005 and these initial values as ^(— 4/xo) = —1 and ^(0) = 0, where the values 



have no particular meaning and are selected to make the plots of the wave function visible. Figs. 1(a) and 1(b) show 
the wave function ^(/i) as a function of /i = 4n/io, and Fig. 1(c) is for the logarithmic scale of |^(/i)|. As Fig. 1(a) 



illustrates, the wave function ^(/i) can be regarded as the sampling of a decaying sinusoidal oscillation at the points 
sufficiently near the origin. However, as seen from Fig. 1(b) ^(/i) drastically changes its behaviour at /i ~ 2400/j,Oj 
so that ^(/i) flips its sign at each step and its absolute value grows up very rapidly for /i > 2400/^0 • We can see from 
Fig. 1(c) that |'I'(m)| grows approximately exponentially for fi > 2400/io- 
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FIG. 1: The wave function *(/i) for (a) < fi/fio < 700 and (b) 2100 < fj./fio < 2450, and (c) the logarithmic scale of its 
absolute value are plotted as a function of /i with A/xg = 0.005, where we choose the initial values as $(— 4/io) = — 1 and 
*(0) = 0. 



This behaviour can be deduced from the following simple argument. We have the recursive relation (|3.4p . For 
/J, ^ ^0) we have 



*(/i + 4mo) - ( 2 - pinl'A ] ^ifi) + ^{fi - 4fio) = 0. 



(3.7) 



If |//^qA/6| ^ 1, it is easily found that the solution is given by 



*(^o+4n^o) ~ An + B, 



(3.8) 
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where A and B are constants. However, if |/i/iQA/6| ^ 1, in the present setting, we generaUy get a solution for which 
the last term on the left hand side of Eq. (I3.7P can be neglected, i.e.. 



*(At + 4^o) 



2 - gMMS^ 
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(3.9) 



For |/i/iQA/6| ^ 1, the factor (2 — /i/ipA/B) is negative and positive if the cosmological constant is positive and 
negative, respectively. In both cases, its absolute value is much larger than unity from the assumption. Thus, ^(/i) 
increases its absolute value at each step by the factor which is much larger than unity and the growth is approximately 
exponential. It is clear that the latter case cannot be regarded as smooth in spite of the very large volume. We can 
expect that this peculiar behaviour becomes prominent for fi ~ 9/iQ ^A~^ for A > and this is consistent with the 
numerical result plotted in Fig. [TJ 

On the other hand, for ij, ^ fiQ, the difference equation (|3.4I) can be approximated by the following WDW equation 
under the assumption that the wave function varies sufficiently slowly (see Appendix |B] for a detailed derivation); 



^ (^/AI*(A*)) + ^m'/^AvI/(^) = 0. (3.10) 

The general solution of the differential equation (|3.10l) for /i > takes the form 

*(/i)=Ai"^ C'lAi (-af n) + C2Bi (^af fi) , (3.11) 

where Ci and C2 are arbitrary constants, Ai and Bi are the Airy functions and we have defined ai := (7r/9)7'^/pjA. 

We present the WDW equations and its general solutions for some operator orderings in Appendix El In the large 

1 

«! //, the above function approaches a decaying sinusoidal curve. Thus, the behaviour of the solution, which is plotted 
in Fig. [H to the difference equation (|3.4p is quite different from that of the solution (|3.11|) to the differential equation 
p.lOp even when the volume of the universe becomes large. In other words, in the equi-area discretisation, the solution 
of the difference equation cannot be regarded as a smooth wave function for large /z in the presence of the cosmological 
constant. We refer to this problem as the large-volume limit problem. 



IV. EQUI- VOLUME DISCRETISATION 

A. Resolution of the large volume limit problem 

In this section, we shall consider the large volume limit in the equi-volume discretisation. We should note that the 
Hamiltonian constraint is written as 
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where the ambiguity in the operator ordering is neglected. Since /z is a function of p, l//z"^ becomes also an operator 
in the quantisation. We put 
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i^ = sin — cos"^ — , EE = sin 
2 2 ' 

We first choose the following operator ordering: 



(I) .„,(!) _„,(!) .,,„(!), 



(4.2) 
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and we will see the other choices in the next section. ^^_^ 

As we have seen in Sec lIIBi in this case, it is convenient to use a new label v instead of /i. We rewrite l/p,^ in 
terms of V as 



1 / 6 



3/2 



KV 



^3 \8ttjIIJ 73 



(4.4) 



by using Eqs. (I2.11[) . (|2.19l) and (|2.24[) . Then, the Hamihonian constraint yields the following difference equation: 

|z; + 4| ||u + 3| - |u + 5||$(u + 4) 



2|w|||w-l|-|i; + l||- 



l^^^%,A\v\]^v) + \v- 



4|||w-5|-|u-3||$(w-4) = 0, 



(4.5) 



where ^{v) — (w|$)- We numerically obtain the solution of the above difference equation with \/3A — 0.1 from the 
initial values $(—4) = —1 and $(0) = 0. Again, we have chosen these values to make the plots visible. The obtained 
wave function <l>(u) is plotted in Figs. 2(a) and 2(b) We see from these figures that the discretised wave function 



oscillates and can be regarded as the sampling of a smooth sinusoidal curve for large v. This means that the large 
volume limit problem in the equi-area discretisation is resolved in the equi-volume discretisation. Indeed, as presented 
in Appendix [B] the difference equation (|4.5p can be approximated by the WDW equation 



^i\-M-)) 



47r 
81K^ 



-fXiM^mv) = 0, 



and the general solution is written as 



^(v) = - (CW^^^ + C2e-*v^") 



(4.6) 



(4.7) 



where Ci and C2 are arbitrary constants and we have defined a2 — of j^a 7'^^piA. The wave function plotted in Fig. [2] 
has good agreement with the solution (14.71) for large v. 



O 





-0.002 



(a) 



(b) 



FIG. 2: The wave function $(u) for (a) < w < 400 and (b) 9400 < v < 10000 are plotted as a function of ?; = 4n with 
V3A = 0.1, where we choose the initial values as $(—4) — —1 and "I'(O) = 0. 



We can also analytically show that the large volume limit should be resolved in this discretisation. For w ^ 1, we 
can approximate Eq. (j4.5p by 

(4.8) 
(4.9) 

(1 - CA) ± J (1 - CA)2 - 1. (4.10) 



$(w + 4) - 2(1 - CA)$(w) + $(w - 4) = 0, 
where C — 2/{27K'^). The characteristic equation for this recursive relation is the following: 



x^ -2{l-CA)x + l =0. 



The roots are given by 



From the above, we can deduce the large-w behaviour of $(w). For < CA < 2, $(4n) purely oscillates sinusoidally. 
For CA < or CA > 2, it is generally dominated by the exponential growth. For CA = 0, $(4n) = An + B, where A 
and B are constants. For CA = 2, $(4n + 4) + $(4n) = (-1)"!?, where D is a constant. Thus, if |CA| < 1, which 
is realistic in our Universe, $(4n) changes very slowly at each step whether it oscillates with the angular velocity 

^ V CA or grows exponentially with the growth rate ^ V CA. This means that we can physically regard $(u) as a 
smooth wave function in the large volume limit. 
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B. Absence of singularity and operator ordering 

In this section, we shall discuss the initial singularity in the equi-volume discretisation. We here demonstrate the 
following four typical types of the operator ordering in the Hamiltonian constraint in the equi-volume discretisation 
and show that the existence or absence of the initial singularity depends on the choice of the operator ordering. 



1 / 96i(sgn(p)) - 1 — , ,_2, 



es, = ,,d^ ^5Sr^^£2+2,^Ai>l. ("1» 



SSi.^^r-^^l^^^rE,2-,'Av). (4,n., 



1 / 96z(sgn(p)) ^l 






We call the above four types of the orderings (a), (b), (c) and (d). Note that EE and l/jj? are commutative with 
each other. We can also take different orderings. For example, l/^i^ could be divided into two and put on both sides 
of F. We will discuss these orderings later. 

The above Hamiltonian constraints for orderings (a) - (d) yield respectively the following difference equations: 



|w + 4|||t; + 3| - \v + b\\^{v + 4) 

L287r 



2\v\ ||« - 1| - Iw + 111 - i^7^-feA|t;| ) $(w) + |w - 4| ||i; - 5| - |t; - 3|| *(w - 4) = 0, (4.12a) 



2|w| ||w + 1| - |w - 1|| - i^73-^A|t;| 1 $(w) + \v\ \\v - 3| - |f - 5|| ^{v - 4) = 0, (4.12b) 



|w| ||w + 5| - |u + 3||$(w + 4) 

1 

"81 ' X2 
|w| ||w + 1| - |w- l||$(u + 4) 

2|«| ||« + 1| - |« - 111 - l^^siei^i^A ^(^) ^ 1^1 11^ + i| _ |„ _ 111 $(„ _ 4) = 0, (4.12c) 

|w + 4|||w + l|-|w-l||$(w + 4) 

2|«| lit; + 1| - |« - 111 - i^^siei^i^A ^(^) + 1^, _ 4| 111; + 1| - |i; - 111 $(« - 4) = 0. (4.12d) 

The WDW equations obtained by the approximation of the difference equations (|4.12ap - (|4.12dp are presented in 
Appendix IB] We can see that $(0) disappears in any of these equations. However, only in Eq. (I4.12al) . $(4) can be 
directly related to $(—4). This leads to the following conclusion. Suppose that we are given $(w) for all w > 0. Then, 
only ordering (a) among the above four choices determines $(u) for all u < 0. For the other choices, the difference 
equation does not determine $(4n) for n = —1, —2, —3, .... In this sense, the absence of the singularity is possible in 
and only in ordering (a) among the choices (a)-(d). Thus, among the six models demonstrated in the present paper, 
the model with the equi-volume discretisation and ordering (a) is the only model where the initial singularity is absent 
and has the large volume limit as a smooth wave function. In ordering (a), if only we fix $(4), we can determine 
$(471) for n = ±1, ±2, • • • . As in the equi-area discretisation, although $(0) is left undetermined, this does not affect 
the unique quantum evolution beyond v — Q. 

Moreover, we can find a general ordering rule required for the absence of the singularity : If the volume operator 
V or its positive power, which is included in EE and l//i^, appears in front of F, the initial singularity appears. We 
can understand this fact as follows. For simplicity, we here assume that there is a constraint equation as C = FV = 
classically. In quantum theory, there are two choices of the operator ordering for the constraint, and then the action 
of the constraint operators Ci — VF and C2 = FV on the state \v) is given by 

Ci\v) = T>F|i-) = -^(K+4|« + 4)-2K|f)+K-4|f-4)), (4.13) 

C2\v) = i?i/|z;)=-^(|i. + 4)-2|i;) + |i;-4)). (4.14) 
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Then, the constraint equations for ^{v) become 

Vv ($(w + 4) - 2$(u) + $(i, - 4)) = 0, (4.15) 

for Ci, and 

K+4*(f + 4) - 2K*(t') + V,^4^v - 4) = 0, (4.16) 

for C2, respectively. In Eq. (I4.15p . there is the same eigenvalue Vy of the volume operator in each term, and then all 
coefhcients vanish for i; = 0. For this reason, in this case, we cannot determine $(—4) from $(4) and $(0), and do 
not have the unique evolution. On the other hand, in Eq. (|4.16|) . the left hand side does not disappear even for v = 0. 
Thus, we can determine the unique solution beyond v = 0, and therefore there is no singularity. 

V. DISCUSSION 

To make the present analysis complete we also need to discuss the matter Hamiltonian Cmattcr- For simplicity, we 
first focus on the equi- volume discretisation. Matter fields can be incorporated by introducing the dependence of the 
wave function on the matter fields such as $(v, (/))■ Then, we can see that the presence of the matter fields does not 
greatly change the discussion. For example, suppose we choose Eq. (I4.12ap as the gravitational part of the difference 
equation. Then, for $(ii, (/>) we obtain the following evolution equation 

\v + A\\\v + 3\-\v + 5\mv + 4, 0) - (^2\v\ \\v-l\-\v + l\\- ^j^^A\v\^ $(„, 0) 

+ \v-A\\\v-5\-\v-3\\ <Piv - 4, 0) = -48if/pi V^C'mattcr$(«, 0), (5.1) 

where we have used the constraint equation Cgiav + C'mattcr — 0. As is explicitly shown in Appendix [Cl for general 
matter fields with an exceptional case, the right hand side of Eq. (|5.1|) is proportional to $(«,(/)) and the coefficient 
vanishes for v — 0. This means that the presence or absence of the singularity does not depend on the matter fields 
in general. It depends on the operator ordering of the gravitational Hamiltonian constraint. This is also true for the 
equi-area discretisation. 

Whether the wave function recovers the smoothness in the large volume limit strongly depends on the choice of 
the discretisation as well as the form of the matter fields. Actually, this problem is already studied well in detail by 
Nelson and Sakellariadou [IllllJ]- In this context, the present paper has shown that if we only include the cosmological 
constant, the large volume limit is problematic in the equi-area discretisation but this problem is resolved in the equi- 
volume discretisation. In reality, the cosmological constant model with each discretisation falls into some particular 
set of the parameter values in their notation and we can find that clearly our result is consistent with theirs. On 
the other hand, we have generated the solution of the finite difference equation for both discretisations and shown 
how problematic the equi-area discretisation is and how nicely the equi-volume discretisation sorts this out. Since 
this behaviour can be explained completely by taking the large volume limit of the discretised Hamiltonian constraint 
equations and the choice of the operator ordering does not affect this limit, it is clear that the large volume limit is 
insensitive to the choice of the operator ordering. 

VI. CONCLUSION 

We have investigated a homogeneous, isotropic and spatially flat universe with the cosmological constant in the 
context of LQC. In particular, we have studied theoretical ambiguities in the quantisation of the Hamiltonian, which 
arise in the operator ordering and the discretisation and are hard to fix within the framework of LQC. We focus on 
two important features of LQC, the absence of the initial singularity and the large volume limit. We have shown 
that the absence of the initial singularity strongly depends on the choice of the operator ordering. Therefore, the 
requirement for the absence of the singularity can potentially make the ordering ambiguity very small. Furthermore, 
we have found a general ordering rule required for the absence of the singularity. On the other hand, the choice of 
the discretisation is crucial when we consider the large volume limit. We have demonstrated two typical choices of 
the discretisation, the equi-area and equi-volume ones. Then, we have shown that in the former case there arises a 
serious problem in the physical interpretation of the wave function as a continuous function in the large volume limit 
and also that this problem is resolved in the latter case. It is clear that we cannot fix these ambiguities only from the 
present results. However, if we can deduce the physical implications of full LQC even qualitatively, we will be able to 
fix the ambiguities of LQC as the phenomenological realisation of full LQG in the minisuperspace model. 
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Appendix A: Wheeler-De Witt equation and its general solution 



Here, we summarise the WDW theory for the flat FRW universe with the cosmological constant. In terms of the 
connection-triad variables c and p, the Poisson bracket takes the form 



{c,p} = 



87rG'7 



and the Hamiltonian constraint is written as 



a 



167r7 



i^(-6c2/H + 2fAb|^ 



0. 



(Al) 



(A2) 



To quantise the system, Dirac's method is used in the WDW theory. That is, the Poisson bracket (IA1[) is replaced 
by the commutation relation between the operators corresponding to the canonical pair: 

8i7r^pj7 



[c,p] 



(A3) 



and the classical constraint (|A2[) is promoted to a constraint for physical quantum states: Cgrav^ = 0. This quantum 
constraint is called the WDW equation. We here consider the following three types of the operator ordering for the 
Hamiltonian constraint: 



ci 



1 



«■■"" 167172/2 



PI 



-6d^^\+2j^A\p\^), 



^^-^ = 16^ {-<^Vmc' + 2^'A\p\^) , 



a(3) = 



-6d^/\p\7;+2YMp\- 



(A4a) 
(A4b) 
(A4c) 



«■•"" 167172/2^ 

If we take the ordinary Schrodinger representation, where p and c respectively act as multiplication and differentiation 

8i7r?|,[7 d^ 



3 dp 
the quantum constraints (jA4a|-(|A4c| yield the corresponding WDW equations 






a(?L*(^) = ^/H 



(7(3) \],(3) „ Stt/^i 



d2^ sgnW d*W _ 1 3 )\ ^^1 ,J^a) ^ 
1^1 dp2 + y^ dp 41^1 "^ j+87r/2^lPl ^ "' 

d2^(2) A 



dp2 



87r/2j 



\p\H^^^^O, 



\P\ ^^2 



^(3) sgn(p)d*(3)\ 



A 



|p|l^(3) =0. 



3 y^' dp2 ■ 2/^ dp J ' Snll^ 
The general solutions of the WDW equations (|A6aP - (|A6cp for p > are given by 



^(1) 



P 



CiAi (-a|p + CsBi (-a|p 



*(2) = CiAi (-alp) + CaBi (-alp) , 



Vl>('^-' — p4 



CiJ, 



8 V3 



aaP^ + C2J 1 -VosP- 



« V3 



(A5) 

(A6a) 
(A6b) 
(A6c) 

(A7a) 
(A7b) 

(A7c) 
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where Ci and C2 are arbitrary constants, Ai and Bi are the Airy functions, J^i is the Bessel function, and we have 
defined as :=3A/(87r?2j)2. 

Appendix B: Large volume limit of the difference equations 



Here, we consider the large volume limit of (/i|Cgrav|*), and then we shall show that it corresponds to the WDW 
equation. First, we consider the operator orderings (jS.Sp and (j3.5|) in the equi-area discretisation. The operator 
ordering p.Sp yields 



(A^ICgravl*) 



1 



16712^2^^2 87r7/2j/i3 



2\V^+^,-V^ 



167r73;2j^3 



^V^ *(m) + |l^M-3/.o - "^M-5/.o I *(M - 4/10 ) 



For ^^ ^Q, we can expand (V^+l^^q ~ ^+M/jo) around /ig as follows: 



-(^)"^^^^{|-<n(^^--^)(^ 



16 V Ai 



O 



(Bl) 



(B2) 



Substituting Eq. ((B2|) into Eq. (jBTj) . we obtain 

1 3 



(Ml^gravl*) 



16^2;2^^2 87r7Z2^^3 \^ 6 



(87n3ny^\,3/2 






M 



-6(^)^v,(^ + 4mo) + ^(a.-4/.o)}+ ''"^3^^'^° AvI.(m) + o[(^^' 



(B3) 



Assuming that the wave function varies sufficiently slowly and expanding the wave function ^(/i±4/io) around /z, we 
obtain 



d*fu) 
vI/(M±4/.o) = *(A*)±-r^(4Mo) 



d/i ''^"' ' 2 dAt2 
Substituting Eq. ((Bit into Eq. ((B3t . we obtain 



i^'*w,,,,,,4i!lw„,^3„of..^'*w 



6 d/x^ 



d/i4 



(M|Cgrav|*> 



n/3 



•^/tt^t^; 



PI 



^(Vm*(/.)) + ^;.3/2A^(^^ ^ o(/.o) 



(B4) 



(B5a) 



Because /i oc p, this equation is just the WDW equation as in Appendix [X] Similarly, we consider the large volume 
limit of (/ijCgravl*) in the operator ordering p.5p . and then we obtain the WDW equation 



(MiCgravl*) 



a/3 



a/tt^t^/ 



PI 



v^|^*(m) + ^t^'/'A^i,.) + 0(mo) 



(B5b) 



Second, we calculate the large volume limit of (w|Cgrav|^) in the operator orderings (|4.12a[) - (I4.12dp in the equi- 
volume discretisation. To do this, we here assume that \v\ 3> 1 and $(w) varies sufficiently slowly. Then, we obtain 
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the WDW equations from the operator orderings (|4.12ap - ()4.12dp . respectively, as follows 

27 



8?pi73/2 V 67r 

27 rs 



8/pi73/2 V 67r 
27 



K 



K 



■^K 



8Zpi73/2 V 67r 

27 rs 



d^ /I 1^/ NN 47r 3 Z?,, ,, . , , ^ /d3$(w) 



, , d^ , , 47r oZ?,, ., , , , /d3$(u) 



8?Pi73/2 V 67r 



X 



ly dw3 

i\vmv)) + ir^'^A\vmv) + o 



d n 1^/ NN 47r . ??,, ,, . , , ^ /d3$(u) 



81 ' K^ 



\ d'y3 



(B6a) 
(B6b) 
(B6c) 
(B6d) 



Appendix C: Matter Hamiltonian constraint 

Here, we shall discuss the matter Hamiltonian operator. We assume that the matter Hamiltonian constraint can be 
written by a'^e{a, </>) for arbitrary matter fields (/>, where r is a constant and e(a, (j)) is a function of the matter fields 
and the scale factor such that e(a, (p) has a nonzero and finite limit for a — > 0. 

First, we consider the case for r < 0. Classically, the matter Hamiltonian diverges for a — >■ due to the inverse 
scale factor. In LQG, such a divergence can be regularised by Thiemann's prescription |18j : we multiply the matter 

Hamiltonian constraint by 1™ = I det e\/ \J\ det £^| 1 , where m can be chosen such that we obtain the positive power 

of the volume factor. Similarly, we regularise the inverse scale factor in LQC as follows J^]. The classical Poisson 
bracket takes the form 



|c,l/3 I =sgn(p) 



%TT-fGl 



\P\ 



l-l 



(CI) 



where we have used V = |p|^^^, and I is the ambiguity parameter similar to m in LQG. Notice that if we choose 
< / < 1, the right hand side denotes the inverse power of p, while the left hand side involves the positive power of 
the volume. Using this property, we rewrite the inverse scale factor a^^ as 



v^- 



3sgn(p) 



{c,vf} 



'\p\ V 87r7G'/ 
where we used \p\ — Vq^^o? . Then, the matter Hamiltonian becomes 

1 



, -, \ 2(1-1) 1 



' 



(C2) 



a2/3„2 



Gr, 



a^e{a, (f>) 
3sgn(p) 



Vo ^e(a, 



{..vf} 



Vo 'e{a,, 



Sn-fGl 
This classical formula can be represented exactly in terms of holonomies as follows: 



C„ 



sgn(p) 
A-K^GlpL 



T.(y:r^h!f{h!f-\V^} 



Vo ^e(a. 



(C3) 



(C4) 



and we can quantise this Hamiltonian immediately by replacing the Poisson bracket {•, •} with ~i[», •]. Using G = /pj 
and Eq. ()4.4p . after some calculation, we obtain 



C„ 



3i(sgn(p)) / 6 



47r7;2j? V87r/pi7 



1/2 



KV 



_ (^sm(-JF3cos(-j-cos(-j^3sm(- 



(C5) 



Then, the matter Hamiltonian operator acts on \v) as 
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Cm&tter\v) — 



,3/2 



3sgn(p) / 6 



SttjII^I 



'pi7 



1/2 



V3 



V- 



V [ 



-211-^ 



Fo ^eo(a, </>)!«) 



(C6) 



where K = (87r7Z|,j/6)'^'^ \'v\/K and eo(a, (/>) is an eigenvalue of e{a,4>). Therefore, \v) is an eigenstate of Cmatter, 
and its eigenvalue vanishes for v = Q. Therefore, substituting this into Eq. (|5.ip . since the eigenvalue of the matter 
Hamiltonian operator vanishes for u = 0, we can uniquely determine the coefficients $(i),(/)) except for $(0,(/)). In 
conclusion, the matter Hamiltonian operator does not affect the absence or presence of the singularity for r < 0. 

Second, we consider the case with r = 0. In this case, the matter Hamiltonian is given by Cmatter = e(a, (j)). Thus, 
the matter Hamiltonian operator acts on the state \v) as 



Cmatter I W) = eo(a, 0)|u). 



(C7) 



Then, \v) is an eigenstate of Cmatter but in this case its eigenvalue does not vanish even for v — Q. Substituting this 
into Eq. (j5.1|) . we find that we cannot uniquely determine ^(O,^) from this equation and hence cannot obtain the 
unique wave function $(u, 0) beyond u < 0. Therefore, in this case, even though there is no initial singularity without 
matter fields, there appears the initial singularity due to the presence of the matter fields. 

Finally, for r > 0, the matter Hamiltonian becomes Cmatter — V^'^ei^a^cf)). Thus, the matter Hamiltonian operator 
acts on the state \v) as 



Cmatter I W) = K^ Eq (fl, 0) | w). 



(C8) 



Therefore, \v) is an eigenstate of Cmatter, and its eigenvalue vanishes for v — Q. Similarly to the case of r < 0, the 
matter Hamiltonian does not affect the absence or presence of the singularity. 
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